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Abstract 

We present a probabilistic proof of the mean field limit and propa¬ 
gation of chaos iV-particle systems in three dimensions with positive 
(Coulomb) or negative (Newton) 1/r potentials scaling like 1/iV and an 
A-dependent cut-off which scales like In particular, for typ¬ 

ical initial data, we show convergence of the empirical distributions to 
solutions of the Vlasov-Poisson system with either repulsive electrical 
or attractive gravitational interactions. 


1 Introduction 

We are interested in a microscopic derivation of the nonrelativistic 
Vlasov-Poisson system. This equation describes a plasma of identical 
charged particles with electrostatic or gravitational interactions 

dt/+p-V,/ + (fc*pt)-Vp/ = 0, (1) 

where k is the (Coulomb) kernel 

cr = {±l}) (2) 

and 

Ptiq) = P[ft]{q) = J fit,q,p)d^p (3) 

is the charge density induced by the distribution f{t,p,q) > 0. 

Units are chosen such that all constants, in particular the mass and 
charge of the particles, are equal to 1. The case cr = -1-1 corresponds to 
electrostatic (repulsive) interactions, while cr = — 1 describes gravita¬ 
tional (attractive) interactions. In the gravitational case, (1-3) is also 
known as the Vlasov-Newton equation. For simplicity, we focus on 
the 3-dimensional case, generalization to arbitrary dimensions d > 2 is 
straightforward. 

* dustin.lazarovici@live.com 
^ pickl@math.lmu.de 
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1.1 Previous results 


While the existence theory of the Vlasov-Poisson dynamics is well un¬ 
derstood - we will cite the pertinent results below - its microscopic 
derivation is still an open problem. To our knowledge, the first paper 
to discuss a mathematically rigorous derivation of Vlasov equations is 
Neunzert and Wick, 1974 m- Better known are the publications of 
Braun and Hepp, 1977 [5] and Dobrushin, 1979 [3], as well as the later 
exposition of Neunzert, 1984 M- For a general overview of the topic, 
we refer the reader to the book of Spohn [TB] . 

Rather than the Vlasov-Poisson equation, the papers of Neunzert, 
Braun and Hepp and Dobrushin treat simplified models with Lipschitz 
continuous forces k G W^’°° = {fc G : ||/c||oo + ||Vfc||oo < oo}. 

The last few years have seen great progress in treating mean field 
limits for singular forces up to but not including the Coulomb case. In 
particular, Hauray and Jabin, 2013, discuss force kernels bounded as 
|^('?)| < with a < d— 1 in d > 3 dimensions [B]. For 1 < a < d— 1, 
they perform the mean field limit for typical initial data and an N- 
dependent cut-off that can be chosen as small as N~^ for a ^ d — 1. 
For a < 1, they are even able to prove molecular chaos with no cut¬ 
off at all. Unfortunately, their method fails precisely at the Coulomb 
threshold a = d — 1. 

In contrast, Kiessling, 2014 proves a non-quantitative approxima¬ 
tion result including the Coulomb singularity under the assumption of 
an (uniform in N) a priori bound on the microscopic forces. The status 
of this assumption, however, whether it is satisfied for generic initial 
data or not, remains open m- 

Recently, Boers and Pickl proposed a novel method for deriving 
mean field equations which is designed for stochastic initial conditions, 
thus aiming directly at a typicality result. With this method, they 
were able to improve the cut-off near the Coulomb case to ~ N~d [1]. 

The aim of this paper is to extend the method of Boers and Pickl to 
include the Coulomb singularity in the large N limit, thus aiming at a 
microscopic derivation of the Vlasov-Poisson dynamics. The Coulomb 
case is qualitatively different from the previously treated interactions 
since the mean field force k*p is no longer Lipschitz, in general, even for 
bounded p. However, we will show how this critical case can be treated 
by exploiting the second order nature of the dynamics and introducing 
an anisotropic scaling of the relevant metric. Moreover, we optimize 
the method in such a way as to achieve a rate of convergence that can 
be faster than any inverse power of TV, depending on decay properties 
of the initial distribution /q. 

An alternative proof, based on similar modifications of the Wasser- 
stein distance, is simultaneously proposed in m- The result presented 
here, however, allows for a significantly smaller cutoff iV“'3+'^ to be 
compared with TV”‘*(‘^+ 2 ) . 
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2 The microscopic model 

Since the Coulomb kernel is strongly singular at the origin, we will 
require a regularization on the microscopic level. For N G N and 
5 > 0, we consider 


f -rV , if |g| > iV-5 

Nf \ I l9l ’ - 


(4) 


On \ {0} this converges to the Coulomb kernel ([5]) as —)> oo. 
Of course, the A^-dependence of the force thus introduced is a tech¬ 
nical necessity rather than a realistic physical model, though similar 
regularizations are commonly used in numerical computations. 

In the mean field scaling, the equations of motion for the regularized 
A^-particle system are given by 


r%(i) =Pi{t) 

I ^ (5) 

]Mt) = Ti E 

for i G 1, Since the vector field is Lipschitz for fixed S,N, we have 

global existence and uniqueness of solutions and hence an A'^-particle 
Hamiltonian flow which we denote by si ^)t ^si^)) ^ 

X . Introducing the A^-particle force — >• given 

by 

1 ^ 

{K^{qi,..,qN))i:=—'^kfiqi-qj), i = l,..,N, (6) 

we can also characterize as the solution of 

^(d/i,(Z),vI/2^(Z)) = (d/i^(Z))), = Z. (7) 

Finally, ii^'^t.oiZ) = (gi(t),Pi(t))i=i,.. _Ar, we define the corresponding 
microseopic or empirical density by 

1 ^ 

H^Z] = nn'^tAz)] := (8) 

Our aim is to show that for typical Z, the empirical density pi^\Z\ 
converges to a solution ft of the Vlasov-Poisson equation as A^ —>■ oo. 


Of course, more general cut-offs can be considered. In the literature, 
the following nomenclature has been established (see e.g. my- 

Definition 2.1. A pair-interaction defined by a kernel A: : R^^ —>■ R'^ 
satisfies a -condition, if 

(5“) 3c > 0,Vg e R'^ \ {0} \Hq)\< |Vfc| < 
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Introducing a cut-off of order N ^ near the origin, the regularized force 
kernel satisfies a {S^)-condition if 

i) k satisfies a (5'“) condition, 

(l) = Hq) for 111 > 

in) \k^{q)\ < N^°‘ for all |g| < N~^. 

In addition, we shall require that 

iv) iVfcf (g)| < for all \q\ < N-\ (9) 

which assures that the regularization around the origin is not too er¬ 
ratic. 

Within this setting, we thus consider 3-dimensional force kernels sat¬ 
isfying a (Sg) condition with a = 2 and the additional assumption 
iv). The lower bound on the cut-off will later be determined as d < f. 
Moreover, we shall adopt the convention fc^(O) = 0, meaning that the 
microscopic dynamics do not contain self-interactions. The reader is 
free to think of ([4]) as defining the microscopic model or consider an¬ 
other regularization of his liking that satisfies the above assumptions. 

2.1 The regularized Vlasov-Poisson equation 

For any <5 > 0 and N G NU{oo}, we also consider the corresponding 
mean field equation 

dtf+p- V J * pt) ■ Vp/ = 0. (10) 

For (formally) N = oo, this reduces to the Vlasov-Poisson equation 
H]). For a fixed initial distribution /o G L°“(]R^ x R^) with /o > 0 
and f f = 1 we denote by fj^ the unique solution of (flUll with initial 
datum /t^(0, •,•) = /q. 


2.2 Method of characteristics 

It is convenient to consider the characteristic flow of the mean field 
system. For N G N,S > 0 and p G L^(R^), we define Kf{--,p) : 
X —>• X R^ by 

K^{q,p\p) ■■= {p,k^ * p{q)). (11) 

Then, the (regularized) Vlasov-Poisson equation (IIOII with initial /o is 
equivalent to the following system of integro-differential equations: 

/o) = /o); Pt) 

pfiQ) = I f^it^Q:P)^^P , . 

( 12 ) 

r(t,.) = <p5(-;/o)#/f 

^‘Ps.si^'Jo) = Z- 


4 




Here, denotes the image-measure of / under ip, defined by 

ip#f{A) = f{ip-^{A)) for any Borel set H C 

In other words, we have non-linear time-evolution in which ip^s {-; /o) 
is the one-particle flow induced by the mean field dynamics with ini¬ 
tial distribution /o, while, in turn, /o is transported with the flow ip^g. 
Due to the semi-group property o (p^, ^ = (p^^ it generally suffices 
to consider the initial time s = 0. 

The method of characteristics can also be though of as establishing a 
kind of duality between the (rescaled) Newtonian dynamics ([S]) and 
the Vlasov equation uni). Indeed, observing that the microscopic force 
can be written as 

i=i 

one easily checks that solves (O with 4'o(^) = 0 if and only if 

gt = fJ.Q is a weak solution of (ITUl) with go = [Z], 

This relation is often used to translate the microscopic dynamics 
into a Vlasov equation, allowing to treat [ZJ and ft on the same foot¬ 
ing. Here, we will go the opposite way, so to speak, and transform the 
mean field dynamics into corresponding JV particle dynamics. To this 
end, we consider the lift of ip^si') to the V-particle phase-space, which 
we denote by That is, for /o G and Z = {qi,Pi)i<i<N, 

we define 

^<i>t,s{Z;fo) ■■= (14) 

Denoting by K : —>■ the lift of the mean field force to the 

V-particle phase-space, i.e. 

(Kt{Z)). ■■= kf * p[f^]{zi), Z = (21,..., zn), (15) 

the flow ^^t,s{Z) = j,(V), g(V)) can also be characterized 

as the solution of the non-autonomous differential equation 

= ($2,(Z),7?,($),(Z))), $,_,(Z) = Z (16) 

to be compared with ©• Finally, we introduce the corresponding 
empirical density 

p^[^t,o{Z)]=ipf^,#p^[Z], (17) 

pertaining to the mean field dynamics with (random) initial conditions 

z e 

In summary, for fixed /g and V £ N, we consider for any initial config¬ 
uration Z £ two different time-evolutions: o(Z’), given by the 

microscopic equations (O and ^^t,o{Z), given by the time-dependent 
mean field force generated by . We are going to show that for typical 
Z, the two time-evolutions or close in an appropriate sense. 
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3 Existence of solutions 


For the well-posedness of the Vlasov-Poisson system, we can rely on 
various results establishing global existence and uniqueness of (weak 
and strong) solutions under fairly mild conditions on the initial con¬ 
figuration /o (Pfaffelmoser, 1990 [111, Schaeffer, 1991 [T7j, Lions and 
Perthame, 1991 [12], Horst, 1993 [Sj). For our purposes, the following 
existence result due to Lions and Perthame is particularly useful: 

Theorem 3.1 (Lions and Perthame). 

Let /o > 0, /o G L^(M^ x R^) fl L°°(]R^ x R^) satisfy 

J\prfo{q,p)dqdp <+00, (18) 

for all m < mo and some mo > 3. 

a) Then, the Vlasov-Poisson system defined by equations (1-3) has 
a continuous, bounded solution f(t,-,-) € C'(R“'"; (R^ x R^)) n 

L°°(R+; L°“(R^ X R^)) for 1 < p < oo satisfying 

sup [\prf{t,q,p)dpdp<-\-oo, (19) 

tG[0,T] J 

for all T < oo,m < mo- 


b) If, in fact, mo > 6 and we assume that fo satisfies 

esssup{/o(g'-f pLp') : \q-q'\ <Rf,\p-p'\ < Rt} 
eL-((0,r)xR3;Li(R3)) 

for all R > 0 and T > 0, then 

sup ||pt(g)||oo < + 00 , VT e (0,-foo). (21) 

tG[0,T] 

Under the assumption of part b) of the theorem, the uniqueness result 
of Loeper, 2006 [Hj then shows that for any T > 0, said / is the unique 
solution in the set of bounded, positive measures on [0, T) x R® sat¬ 
isfying /|j_g= /o in the sense of distributions. Moreover, it is known 
that as long as the charge density is bounded, solutions with smooth 
initial data remain smooth (see e.g. in [7|). 

As Lions and Perthame remark - and as one can verify by following 
their proof - part b) of the theorem actually yields a bound on the 
charge density that is uniform in N if one considers a sequence of 
regularized time-evolutions as (for instance) in dTIIll . We will note this 
important fact in the following lemma. 

Lemma 3.2. Let fo £ L^(R^ x R^) fl L°°(R^ x R^) and f^ he the 
solution of the regularized Vlasov-Poisson equation (HOI) (with corre¬ 
sponding cut-off) and initial datum /^(O, •,•) = /q. If fo satisfies 
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assumption (EOD of the above theorem, there exists a constant Cp > 0 
such that 

llpflloo + llpnii < Cp, VA^eNU{oo}, Vt>0, (22) 

where and, with a little abuse of notaiton, pf^ = p[ft]- 

Since condition (1201) is rather abstract, we want to state a more intu¬ 
itive sufficient criterion. 

Lemma 3.3. Let fo G x K^) (~l x R^), / > 0. Sup¬ 
pose there exist functions p € L°°(R^') and i^dpl) S with § 

monotonously decreasing and an S > 0 such that for all |j>| > S 

fo{q,p) < piq)'i^i\p\)- 

Then /g satisfies assumption (Enn- Special cases: 

• /o has compact support in the p-variables. 

• fQ is a thermal state of the form p{q) e~^^ with\\p\\oo< oo, fi > Q. 
Proof. For given R,t > 0 we have to consider the function 

f{t,q,p) ■■= esssup{fo{q -\-pt,p ) : \q - q\ < \p-p \ < Rt}. 
Choosing R' > S RT, we have 

J f{t,q,p)d^p= J + J 

<|7ri?'^||/(t,-,-)lloo + IIpIIoo /" sup d{\p'\)d^p 
u J |p— 

<^7ri?'3||/o||oo + IIpIIoo j ^{\p\-Rt) d^p 

^Cll/olloo + ||p||oo|ll?||l < OO, 

where in the second to last line we used the monotonicity of i?(|p|) and 
the fact that ||/||oo = ||/o||oo- □ 


One crucial consequence of the bounded density is that the mean field 
force remains bounded, as well. 

Lemma 3.4. Let k be the Coulomb kernel, and p G L^ r\L°°{ 

Then there exists C > 0 such that 

||fc*p||oo<C||p||}/^||p||L/^ 


(23) 


Proof For R > 0, we compute: 
||fc*p||oo< J k{y)p{x-y)d^y 


Ky)p{x - y) 


\v\<R 


\v\>R 


< 


J + ^IIpIIi =47ri?||p||oo+ i?' 




This last expression is optimized by setting R = (47r) || p|| 

which yields ||fc*p||oo < 2 ( 47 r) 2 C||p||C3||^||^3^ 




□ 
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4 Statement of the results 


In the following, all probabilities and expectation values are meant 
with respect to the product measure given at a certain time by 
That is, for any random variable H : —>■ R and any element A of 

the Borel algebra 


Ar 



(24) 



(25) 


Note that since ^^t,s leaves the measure invariant 





In particular: 


(Z e A) =F^G A). 


(26) 


We will often omit the index N for Pg and Eg defined with respect to 
the product measure 

To quantify the convergence of probability measures, we will use the 
Wasserstein distances (also known as Monge-Kantorivich-Rubinstein 
distances). In the context of kinetic equations, they were first intro¬ 
duced by Dobrushin in [2. We shall briefly recall the definition and 
some basic properties. For further details, we refer the reader to the 
book of Villani |201 Ch. 6]. 

Definition 4.1. Let P(R"’) be the set of probability measures on 
R". For given G P(R"), let n(/r, i^) be the set of all probabil¬ 
ity measures R” x R" with marginal /i and i/, respectively. Then, for 
p G [1, +oo), the p’th Wasserstein distance on P(R”) is defined by 



(27) 


B^xR" 


Convergence in Wasserstein distance implies, in particular, weak con¬ 
vergence in 7^(R"), i.e. 



k —>■ oo, 
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for all bonded, continuous functions $. Moreover, convergence in Wp 
implies convergence of the first p moments. Wp satisfies all properties 
of a metric on P(R"’), except that it may take the value +c». 

The most common version is the first Wasserstein distance, for 
which we have the Kantorovich-Rubinstein duality: 

Wi{p,iy)= sup \ [ g{x)dp{x) - [ g{x)di^{x)], (28) 

where ||g||Lip := sup for g : R" —K. We will also consider 

x,v 

the infinite Wasserstein distance defined by 


Woo(Aij ^) = inf{7r — esssup |a; — j/| : tt S n(/r, z/)}. (29) 

We can now state our precise results in the following theorem. 

Theorem 4.2 (Molecular chaos). Let fo G x R^) a probabil¬ 

ity measure satisfying the assumptions of Theorem \3.1\ a) and b) and 
f^ the unique solution of the regularized Vlasov-Poisson equation 
with initial datum /q. For 0 < 5 ^ let be the N-particle flow 

solving with cut-off width N ^ and let be the N-particle 

mean field flow induced by f^ as defined in d. Then, for any 

T > 0, there exists a constant Cq depending on sup ||/9[/^]||/,oo(ro TlxR^) 

atgn 

such that for any /3 > 0 there exists a constant Cp such that for all 
N>No:= 


3t e [0,T] : rd>tp[Z) - ^$t,o(^)|oo > N 


T-S 


< 


TCp 
NP ’ 


(30) 


where j-loo denotes the maximum-norm on 

This result implies molecular chaos in the following sense: 

Corollary 4.3. Let := fo and F^^ := ^d/tpA^Fo the N-particle 
distribution evolving with the microscopic flow O- Then the k-particle 
marginal 


WF,'^(2i,...,Zfc) := J F^^{Z)d^Zk+i...d^ZN (31) 

converges weakly to ft as N ^ oo for all k G N, where ft is the 
unique solution of the Vlasov-Poisson equation O with f^\t=o= fo- 
More precisely, under the assumptions of the previous theorem, we get 
a constant C > 0 such that for all N > No 

Wi{‘'^^Ft^,(^'^ft) < VO < f < T. (32) 

Proof. For fixed 0 < t < T, let A C R®'^ be the set defined hy Z G 
A \^d!tfi{Z) — ^^t,o{Z)\^ < N~^. Hence, according to the 

previous theorem, Po(Vl°) < for sufficiently large N. In view of 
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the Kantorovich-Rubinstein duality (I^Rll . we have: 




= sup 

II^IIltp—1 

= sup 

\\g\\Liv^i 

= sup 

WgWhip^i 


Zk)d^Zi...d^Zk 

F^{Z) - ft{Z))g{zi ,..., zu)d^zi...d^zu...d^ZN 
J {d>t,o#F^{Z) - <!>t,o#F^{Z))g{z,,..., Zk) d^^Z 


Introducing the projection Pk : —t R^, (zi,..., zat) i-)- (zi,...,Zfe), 

this can be further rewritten as 




= sup 

\\g\\Lip^i 

= sup 

\\g\\Lip^i 


sup 

l\\Lip — ^ 


Fo^{Z){giPkd>t,o{Z)) - 9iPkd>tAZ))) d^^Z 
Fo^(Z)(g(Pkd^t,o(Z)) - g(Pkd>t,o(Z))) d<^^Z 
Fo^(Z)(g(Pkd^t,o(Z)) - g(Pkd>t,o(Z))) d<^^Z 


(33) 

(34) 


Using that all test-functions are Lipschitz with UgllLip = 1, we have 
m < Po(^")||i"o'^||oor4/i,o(^)-'^<i>t,o(^)|, with IlFo^lU = (ll/olloo)'^. 
Recalling that |^4'o_o(^) — ^‘bo,o(^)|oo = 0, we have 

t 

\^dil^{Z)-^‘^l^[Z)U < J |iff(4'i,o(Z))-7?($,%(Z))Uds, 

0 

t 

l^d^loiZ) - ^d>loiZ)\oo < J r^^,o(^) - o(Z)Uds. 

0 


The mean field force K is of order 1 ILemma 13.4|) . while the micro¬ 
scopic force is bounded by . Hence, there exists a constant 
C > 0 such that j^'I'j q(Z) — ^4>(_q(Z)|oo < TCN"^^ and consequently 
l^diloiZ) - ^$)o(Z)joo < T^C'N^^ for all t < T. Choosing (3 := 35 
in dsni we thus get another constant C” such that 

(IMl) < C” max{r2, T^}N-\ VO < t < T. (35) 

On the other hand, for Z £ A, we have for any g with ||g||Lip = 1, 

\g{Pu^d>t^Q(Z)) - g{Pk^<^t^^{Z))\ < r'bt,o(^) - '^$t,o(Z)U < iV"^ 

for all t < T and thus (IMl) < N ^. Together with (1551) . we get a 
constant C" such that 

IUi(('=)Ff, oVD < C"'{1 + T^)N-\ VO < t < T. (36) 
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Finally, we will prove in Proposition 19.II that 


ft) < Vt < T (37) 

where Cq depends only on /o and T. Putting everything together and 
using the 

statement follows. □ 

It is a classical result in probability theory (see e.g. [9], [5], |191 Prop.2.2]) 
that molecular chaos in the sense of the previous corollary implies con¬ 
vergence in law of the empirical distribution nflZ] := to 

the constant variable ft- However, under additional assumptions on 
the decay of /o, we can obtain the following quantitative result: 

Theorem 4.4 (Particle approximation of the Vlasov-Poisson system). 
Let /o S X R^) a probability measure satisfying the assumptions 

of Theorem \8.1\ a) and b). For 0 < <5 < ^, let be the N-particle 

flow solving with cut-off width N Let p £ [l,c») and assume 
that, in addition, there exists k > 2p such that Jjjg|z|^d/o(2:) < -|-oo. 
Then, the empirical density ix^[Z] := Hq [i/t,o{Z)] converges to the so¬ 
lution of the Vlasov-Poisson equation in the following sense: 

For any T > 0 and 7 < min{i,^,5}, there exists a constant Cq 

depending on /g and T and constants c,Ci,C 2 depending on k,p,^ 

f 2(CnT+l) Y 

such that for all N > Ni := eV J 


e [0,r] : Wpifx^'lZjJt) > 


(38) 


where f is the unique solution of the Vlasov-Poisson system on [0, T\ 
with /|(^o= /o md 6 V 2p := max{ 6 , 2p}. 

The proofs of the theorems will be given in Sections [5] and [HI 

Remarks 4.5. 

1) Our results allow to choose the width of the cut-off arbitrary close 

to which corresponds to the scale of the typical distance 

between a particle and its nearest neighbor. 

2) The results can be straightforwardly generalized to include ex¬ 
ternal forces or noise, provided that the regularity of solutions to 
the corresponding mean field equation remains as assumed. 

3) All dynamical estimates required for Theorem 14.41 would hold 
directly for the infinite Wasserstein distance. However, there are 
no concentration estimates available in terms of Woo to establish a 
good approximation of the initial /o by the empirical measure p,Q . 
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5 A new measure of chaos 


The strategy of proof is to control the deviation of the microscopic time 
evolution from the mean field time evolution in terms of the following 
-/V-dependent quantity: 

Definition 5.1. Let ^^t,o the mean field flow defined in (ED and 
0 the microscopic flow defined in We denote by = 

{qi{t))i<i<N and Q = {pi{t))i<i<N the projection onto the spatial, 
respectively the momentum coordinates. 


Consider the quantity ^ defined as 


where \Z\ao = max{|zi| : 1 < * < 3A} is the maximum-norm on . 
In view of Theorem 14.21 our aim is to show that for any e > 0: 


sup {Af,J>A 
Lo<s<r 


-(5 


—^ 0 , 


(39) 


faster than any inverse power of N. This will be done by introducing 
for A > 0 and N G N the stochastic process J^’^{Z) given by 

:=min|l, sup |eVi°sW(r-«) (A-^Af, -h 

L 0<s<£ J 

(40) 

and controlling the evolution of (J^ ’ ). 


The crucial innovation with respect to [T] is that distances in spatial 
and momentum coordinates are weighted differently by a A-dependent 
factor (here: y^log(A)), exploiting the second-order nature of the dy¬ 
namics. Moreover, the quantity has been redefined in order to 

optimize the rate of convergence. 


The relevance of (I39II for the proof of Theorem 14.41 is grounded in the 
following observation: 

N 

Lemma 5.2. For X = (xi,...,XAr) G let Pq[X] '■= G 

i—1 

Then we have for all p G [l,oo]; 

Wp{p^[X],p^[Y])<\X-Y\. 


This implies, in particular, for any ^ > 0 


sup Wp{p‘^[^,,o{Z)],pS[<^>s,o{Z)]) > e 

'-0<s<t 


<Po 


sup {Af } > ^ 

L0<s<T 


( 41 ) 
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Proof of the Lemma. Since Wp < Wq for p < g, it suffices to consider 
the infinite Wasserstein distance defined in (l2^ . 

N 

We then observe that ttq = X) G n(p^[Z],p^[F]) with ttq — 


i—1 

esssup \x — y\ = max \xi — yA = \X — Floo- CH 

l<i<Af 

In total, we will split our approximation result for the Vlasov-Poisson 
equation into 

[Z],/t) < Wp{^^^['^>tAZ)],^^o{<^tAZ)]) (42) 

+ Wpifi^[^tAZ)]Jt^) (43) 

+ Wpif,^Jt). (44) 


The first term dUl) is the interesting one, concerning the difference 
between microscopic time-evolution and mean field time-evolution. It 
will be controlled in terms of Eo( and shown to converge in prob¬ 
ability faster than any inverse power of N. 

The second term Wp(At^[$t,o(^)],/t) = ITp(pffo#p^[Z], p^o#/o) con¬ 
cerns the sampling of the mean field dynamics by discrete particle tra¬ 
jectories. We will use a recent large deviation estimate of Fournier 
and Guillin [3] to determine the typical rates of convergence for the 
initial distribution and then control the growth of (H51) by a Gronwall 
estimate. 

Convergence of 631) is a purely deterministic result: solutions of the 
regularized Vlasov-Poisson equation approximate solutions of the 
proper Vlasov-Poisson equation © as the width of the cut-off goes to 
zero. 

The central idea of our strategy is thus to hrst sample the (regularized) 
mean field dynamics along trajectories with random initial conditions, 
i.e. approximate f^^ by /i^[$i o(Z)], and then control the difference 
between mean field trajectories and the “true” microscopic trajecto¬ 
ries in terms of the expectation This approach has several 

important virtues: 

1. The method is designed for stochastic initial conditions, thus al¬ 
lowing for law-of-large number estimates that turn out to be very 
powerful. (Note that the particles evolving with the mean field 
flow remain statistically independent at all times.) 

2. The metric A'^tAZ) ~ ^4>t_o(-^)|oo is much stronger than usual 
weak distances between probability measures, thus allowing for 
better stability estimates. 

3. Since < 0 if sup A'HsAZ) - ^^sAZ)\oo > N-^ we 

0<s<t 

only have to consider situations in which mean field trajectories 
and microscopic trajectories are still close together. 

4. Exploiting the second-order nature of the dynamics, we weigh 
distances in x-space and momentum space differently, with an 
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-/V-dependent factor Y^log(-/V). As we compare microscopic tra¬ 
jectories to characteristic curves of the mean field equation, the 
growth the spatial distance is trivially bounded by the difference 
of the respective momenta. The idea is thus to be a little more 
strict on deviations in space, so to speak, and use this to obtain 
better control on fluctuations of the force. 


6 Local Lipschitz bound 

If all forces were Lipschitz continuous with a Lipschitz constant L inde¬ 
pendent of A^, we could easily conclude that — < 

(1 -I- L)\^'i’tfi{Z) — ^<^tfi{Z)\oo, from which the desired convergence 
readily follows. However, the forces considered here become singular in 
the limit N ^ oo and hence do not satisfy a uniform Lipschitz bound. 
Nevertheless, we observe that, for the mean field force , the 

global Lipschitz constant \\k^ * \ \Lip diverges only logarithmically as 

the cut-off is lifted with increasing N. Due to the pre-factor -^log (N) 
in Definition 15.11 the particular anisotropic scaling of our metric will 
allow us to “trade” part of this divergence for a tighter control on spa¬ 
tial fluctuations. This will suffice to establish the desired convergence, 

using the fact that = N grows slower than N'^ for any 

e > 0. (C.f. also [Tl] where we have implemented the same idea). 


We summarize our first observation in the following Lemma. 

Lemma 6.1. Let 0 < 7 < 1 and assume that Z : — >■ satisfies 

1 ^( 9 )! < c • min{Af^'’', (45) 

for some c > 0. Then there exists a constant Ci > 0 such that 

\\l * Pt(a;)||oo < C/maxIl, \/log(Af)} (||pt||i + ll/OtlU)- (46) 
Proof. We estimate 

IK*Pt(x)||oo = I l{x-y)pt{y)d^y 

J 00 

< / l{x-y)pt{]j)d?y 

J 00 

\x — y\<N~'^ 

+ J Kx-y)pt{y)d^y 

N-^ <\x-y\<l 

+ J l{x-y)pt{y)d^y 

\x-y\>l 

The first term is bounded by 


lix-y)ptiy)d^y < ||pt||ooAr^'^|B(Af ^)| < xTt ||pt||oo, 

00 «J 


\x-y\<N-^ 
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where B(r) denotes the ball with radius r. The last term is bounded 

by 


\x-y\>l 


y)pt{y) 


< c\\pt\\i. 


Finally, the second term yields 


iV-'>'<|x-y|<l 


y)pt(y)d^y < llptll 


Af-T'<|y|<l 



< 47rc7||pt||oolog(iV). 


□ 

One immediate application of the Lemma is to l{q) = {q), showing 

that a regularized mean field force is Lipschitz continuous with a con¬ 
stant proportional to \og{N). Our aim is now to prove that for typical 
initial conditions, fluctuations in the microscopic forces can be bound 
in a similar fashion, as long as and are close. 

Definition 6.2. Let 




^ ,if|g|>3fV-^ 

, else 


(47) 


and L : —)> be defined by {L{Z))i ■= ^ ~ 1j)- 

Furthermore, we define Lt{Z) by {Lt{Z))i := if * = J *q 

f{t,qi,p) dV 

Lemma 6.3. For any ^ G with |^|oo < 2N~^, it holds that 

\kf{q)-kf{q + 0\^<lf{qm^. (48) 

Proof. First note that by assumption the derivative of k^ is bounded 
by so that (H51) holds for |g| < SN~^. For |g| > there exists 

s G [0,1] such that \kf{q) - kf {q + f)\ < |Vfcf (g -f s^)|oo|C|oo, where 

|Vfcf(g + sOU <2|g + sCr"- (49) 


The expression on the right-hand-side takes its greatest value if f is 
antiparallel to q and s = 1. Hence, we have 

\kf{q) - kf{q + OU < 2 |g(l - ||)|-" ICU- (50) 

Since |g| > 3N~^ and |^| < 2N~^, it follows that ||| < §■ Thus, we 
get \kfiq) - kf{q + Oloo < 2(^)' l^loo < ^ l^loo. □ 
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7 Law of large numbers 

In order to control the evolution of Eo(J(^), we will require as an 
intermediate step that the mean field force (and its derivative) can be 
approximated by the analogous expressions for the discrete measure 
with random Z. The key observation here is that if the 
A^-particle configuration evolves with the mean field flow the 

particles remain statistically independent for all t, thus giving rise to 
a law-of-large-numbers estimate. 

Definition 7.1. For any t > 0 and fixed i5 < h, we consider the {N 
and t dependent) sets At,Bt,Ct defined by 

ZGAt ^ \J!''\z)\ < 1 

ZGBt ^ |iFi^($t.o(Z))-7?($i.o(^))L <iV-i+2^ 

Z€Ct ^ |Lf($i,o(Z))-!(<!>*,o(^))L < 1 


where K is the mean field force (TTSl) and L as in Definition 16.21 

We now want to show that for any t, initial conditions in Bt D Ct are 
typical with respect to the product measure Fq := on 

Proposition 7.2. Let pt G L^(R^) fl with ||pt||i = 1 as 

before. Let h : —>■ K and suppose that for given d > 0 and 

N G N there exists c > 0 and an exponent 2 < a < 3 such that 
|/i(a:)| < c • min{A'^“'^, | 9 |~“}, V(? G R®. Assume furthermore that 


6 < min 


in I 


1-2^ 1-/3 


}■ 


.2a — 3’ a 

Then there exists for all j > 0 a constant > 0 such that 
. N 


sup 

‘-l< 2 <Ar 






> N-f^ 


Proof. Let 


N 




N N 

and D := U A- Then ¥{D) < X) P(A) = 3VP(Di). 

i=l z=l 

By Markov’s inequality, we have for every M G N: 

r 1 ^ 

Pt(A) <Et Kqi - Qj) - h * Pt{qi) 


N 


2M-\ 




ij 3 ^E[(^(/l(gi -qj)-h* pt{q^)fj 


2M 




(51) 


< (52) 


(53) 


(54) 
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Let Ai '■= {k G Nq I |k| = 2M} the set of multiindices k = {ki,kN) 

N 

with ^ kj = 2M. Let 


N 

G*" := Y[{h{qj - qi) -h*pt{qi))^A (55) 

Then: 

N 

j=l ksAI ^ ^ 

Now we note that Et{G^) = 0 whenever there exists st 1 < j < N such 
that kj = 1. This can be seen by integrating the j’th variable first. 

For the remaining terms, we have for any 1 < m < M: 

J \Hqi - Qj)rMqj,Pj)d^Pjd^Pj = J \hr{qi - qj)ptiqj)d^qj- 

Now for 2 < a < 3 and m = 1 we estimate 

\h*ptiqi)\< J\h\{qi-y)pt{y)d^y 

<c y \y\~°'pt{qi-y)d^y + c J \qj\~°‘Pt{qi - v) d^y 
|y|<i l!/l>i 

< c(47r||pt||oo + ||pt||i)> 

while for a = 3, we find: 

\h*pt{qi)\< yi^Kgi-2/)pt(2/)d^2/ 

< c( y + y + y 'j\hiy)\ptiqi-y)d^y 

|y|<JV-5 Ar-5<|y|<l |y|>l 

<c||pt||oo y d^y + c\\pt\\oo J j^d^y + c J pt{qi-y)d^y 

|y|<JV-5 Af-5<|y|<l |y|>l 

< C (47r||pt||oo(^ + log(fV'^)) + llptlli) . 

For m > 2, we hnd in any case 

y \hriqi - y)pt{y) d^y = J \hriy)ptiqi - y) d^y 

< y \hr{y)pt{qi-y)d^y+ J \hr{y)ptiqi - y)d^y 

|y|<JV-5 |y|>W-5 

<c||/9t||oo(47riV-3^7V“^™+ y ^d3y)<87rc||pt|UiV(“’”-'^'. 

|y|>JV-5 
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Hence, setting Ca ■= 167rc||/9t||oo (l + l{a= 3 } log(-/V)) we can conclude 
that for all TO > 2 : 

\h{qj - gO - h(g,)r < (57) 

Now, for k = (fci, ^ 2 ,..., fcAf) G A4, let denote the number of kj 
with akj 7 ^ 0. Note that if > M, we must have kj = 1 for at least 
one 1 < j < N, so that Et(G'^) = 0. For the other multiindices, we get 
(using that the particles are statistically independent): 


E 


:(G‘") = Et [^]^(fe 5 (gj -qi)-k* ptiqi))’' 
j=i 
N 

<n]Et {\h{qj - qi)\ + \h * pt{qi)\)’' 


(58) 


(afej—3)5 


N 

t=i 

^^2Af j\^2AIckS 

Finally, we observe that for any ^ > 1, the number of multiindices 
k G with ^k = I is bounded by 

1 < ( 7 ) ( 2 ^)* ^ {2Mf^NK 

#k=i ^ ^ 


Thus: 


HDi) < 


1 


_/y2M(l-/3) 


-/3) X! 


kGVVI 


2M 


Et(G*^) 


< CI^Cm 


N' 


2MaS ^ 


^2M(l-/3) 


-/3) 


-3S)l 


1=1 


< CI^MCmN-^^, 

where Cm is some constant depending on M and 


e := 


1 - 2/3 + <5(3 - 2a) if 3^ < 1 
2(1-/3-0(5) if 3(5 >1. 


(59) 


e > 0 according to m- For 2 < a < 3 we conclude the proof by 
noting that 


Ft{D) < NPt{Di) < C^^MCm 


(60) 


and choosing M so large that {eM — 1) = 7 . For a = 3, however, 
equation (I5H becomes 


t{D) < G'(M)(1 + \og{N)f^ 


(61) 
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where C'{M) is some constant depending on M and ||pt||oo- This can 
be rewritten as 

(1 + \og{N)f^ Af+b (62) 

The function (/(x) = ^ ^ S [1, oo) is continuous with lim g{x) = 

0. Hence, it has a maximum C < +oo. In particular, < C in¬ 

dependent of N and the announced result holds for a = 3, as well. □ 

Corollary 7.3. Let Bt,Ct as in Definition \7.1\ Then we find for any 
7 > 0 a constant such that 

r^iBt) > 1 - 
<(« > 1 - 

In other words, for any fixed t, initial conditions in Bt H C* are typical 
with the measure of “bad” initial conditions decreasing faster than any 
inverse power of N. 


Proof. Note that 


\K^{Z)-KiZ)\^<N-^+^^ 

N 


max 

l<i<N 


^ ^ (ft - ft) - * pf (ft) 




< N 


-1+2(5 


and similarly 


\LfiZ)-LiZ)\^<l 

N 


max 

l<i<N 


Jf - ft) -Is * pf (ft) 




< 1 . 


Applying the previous result once for with a = 2 and /3 = 1 — 2S 
and once for l^ with a = 3 and ,5 = 0, we get 


max 

Ll<i<Af 


max 

Ll<i<Af 


N 


N 


N 


ksA - ft ) -kf * Pt{qi) 


j¥=i 


N 


> N 


-1+2(5 


< 


iV7’ 


A - Qj) - *pf(ft) 


j¥=i 


> 1 


< 


Cj_ 

N'T 


(63) 

(64) 


Observing that = Pf'($t.o(-Bt)) and P^(Ct) = Pf ($t,o(Ct)), 

the statement follows. □ 
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8 Proof of Theorem 14.2 


We now have everything in place to prove our first Theorem 14.21 That 
is, we will prove the following statement: 


Under the assumptions of Theorem \4.<l\ we find for all 6 < ^ and T > 0 
a constant Cq > 0 such that for any j3 > Q there exists a constant C > 0 

, CoT+1 ,2 

such that for N > e^ > 


sup {A 

0<s<T 


}>N 


-s 


< T 


C 

NP' 


(65) 


In order to control the evolution of sup {A^g}, respectively 

0<s<t 

defined in (I4()L we will need the following Lemma. 

Lemma 8.1. For a function g : ffi. —>■ R, we denote by 


d+g{t) := 


lim 

A^t\0 


g{t + A^t) - g{t) 
A^t 


( 66 ) 


the right-derivative of f with respect to t. Let g € ^^(R.) and h(t) := 
sup g{s). Then h{f) exists and h(t) < max{0,g'(t)} for all t. 

Q<s<t 

Proof. We have to distinguish 3 cases. 

1) If g'{t) < 0, there exists A^t > 0 such that g{s) < g{t),ys G 
[t,t + A^t). Thus for all t' G [t,t + A^t) we have h{t') := sup g{s) = 

0<s<t' 

sup g{s) = h{t) and h{t) = 0. 

0<s<t 

2) If g(t) < h{t)^ there exists A^t > 0 such that g{s) < h(t)ys G 
(t— A^t,t + A^t). This means that h is constant on {t—A^t,t + A^t) 
so that, in particular, h{t) = 0. 

3) If g(t) = h(t) and g'{t) > 0, there exists A^t > 0 such that g is 
monotonously increasing on {t — A^t, t-\-A^t). Hence, we have h{t') = 

sup g{s) = g{t') for all t' G [< + A^t) and thus h{t) = g'{t). □ 

0<s<t' 

Proof of Theorem 14.21 Recall from Definition 15.11 

Al, = 

As announced, we introduce for any A > 0 and A G N the process 
Jt'^{Z) := min 11, sup || 

L 0<s<t ^ J 


We consider the expectation which we split as follows: 

Eo(jf’^) = Eo(jfI Al?)+Eo(jf’^ I At\BtnCt)+Eo{Jt^’^ \ AtnEtOCt). 
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Here, At,Bt,Ct are the sets defined in Definition 17.11 and \ A 


denotes the restriction of to the set At C 


etc. 


1) On At, we have = 0, since Jt'^{Z) is already maximal 

and thus also 

|^^)=0. (67) 

2) For Z G At, we have to consider J^'^{Z) < max|o, 
with 




= - Av'log(A^)e^Vi°s(")('r-‘) {N^A^t + 


( 68 ) 




We split dtA^t into 




a*rvi,i (Z) - (^)loo < Idt^Koiz) - ^^Uz))\ 


,0\ 

N ,T, 2 


< r^ioiz) - ^$^.o(^)u < sup r^uz) - ^$f.o(^)ic 

0<s<t 


(69) 


and 


dtr^fioiz) - ^^Koiz)\^ < \dtrKoiz) - ‘^Koiz))\^ 

< |iFf(vI/io(Z))-7?*($io(Z))|, 


(70) 


We begin by controlling the contribution of “bad” initial conditions not 
contained in Bt and Ct- Since \k^\ < , the total force acting on 

each particle (with the 1/iV-scaling) is also bounded as \K^{Z)\ao < 

. The mean field force K is of order 1, according to Lemma 13.41 
and q(Z) — g(Z)|oo < 1 since Z S At- In total, we thus 

have sup{|a+jf’^(Z)| : Z G At} < for some C > 0. 

According to Prop. 17.21 the probability for Z ^ BtHCt decreases faster 
than any power of N. Hence, we can find for any At > 0 a constant 
(which may differ from the in Proposition [721) , such that 


dtW^'^\At\{Btr\Ct)) 

< sup{|9+jf’\Z)| : Z G A}Po[(A nSO" 


< e 


A^/iog(jv)r 


(71) 


iV« 


3) It remains to control the change of for the typical initial con¬ 
ditions Z G At C\BtC\Ct. To this end, we consider: 


\K^{^l,{Z))-Kt{^lt,{Z))U 

< |A'f(Tio(Z))-A'f($ig(Z))|, 

+ \K^{<^lt,{Z))-Kt{<^>i,{Z))\.o 


(72) 

(73) 
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Since Z G Bt, it follows that 

\K^{^loiZ)) - KtKo{Z))\^ < N^^-\ (74) 

For (17^ . we use the triangle to get for any 1 < i < 

{Kf{^l^{Z)) - K{<^l^{Z)))^ < fcf (4/1 - 4,1) _ fcf (4>i - $1) 

i-1 

N 

< - 'I']) - - ‘J’DIoo- 

i=i 

Since Z G At, i-e. J^’^{Z) < 1, we have in particular for iV > 3 i.e. 
log(-/V) > 1, sup |1^4/1 q(Z) — ^4)1q(Z)|oo < N~^. Therefore we can 

0 <B <t 

use Lemma 1(1.31 to get the bound 

|Af (4/1 - 4/1) - fcf ($1 - 4>i)|^ < If ($1 - 4>i)|(4/i - 4/1) - (4/1 - 4/i)U 

<2lf(4>i-$1)14/1 o-$ioU. 

Since Z G Ct, it follows with Lemma EH] that 

N 

E - <&]) = (L^i^tAZ)), < IKf * pf (9)||oo + 1 < 2Q log(7V)(l + llpf lU). 

i=i 

Hence, setting Cq := 2(7; Cp, where Cp is the uniform bound on the 
charge densities from (l22|l . we have 


|-|4/2(Z) - $?o(z)U < Colog(iV) \^IAZ) - ^\AZ)\^ + 


dt 

for Z G At Together with this yields: 

dtKt 


(75) 


AtnBtnc, 


< AA^iAZHoiz) - AAZ)A + ^,Hoiz) - <^lo{zA 


<AAi(N)\^l,{z)-A,AZ)A 

+ [Colog(iV)|4/io(Z) - $io(^)|oo + 

< Co 

Plugging this into equation (1681) and using that C; > d/r we have found 
jf’'^(Z) < max{0, if’"^(Z)} with 

/f’^(Z) < -Av'log(Ai)e^V>°g(^)(^-*) t +7V3 ^-i) 

+ gAViog(Ar)(r-t)^5 (^Cox/IogTVAf t + iV^-^-i 

= A^ogiN)NA^A°s{N)iT-t) 


(Co -A)Af,+ 


1 




-A A^'^-i 
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Choosing A = Cq gives that this is negative. Hence, we have 


d+Eo{J^’^ \AtnBtnCt) = o. 


Together with (1^71) and (I7ip we have found: 


a+Eo(jf’^) < 


( 76 ) 


and thus 

Eo(jf’^) - Eo(Jo'^’^) = Eo(jf- J^'^) < (77) 

uniform in t £ [OiT]. Observing that A^o = 0, we have = 

eW^°s{N)T]^3S-i _ demand 

A> 1. (78) 


(Here we exploited the fact that grows slower than any power 

of N). The random variable jf 

is then certainly non-negative and it follows from (1771) that: 


4’\Z) - J^’\Z) > i' 


< oygA^/loa(Af)T 

- 7V7' 


However, if J^'^{Z) — Jq'^{Z) < i, we have together with (ITRl) that 
J^’^{Z) < 1. And in this case, we can conclude from 

4'\Z) - J^’\Z) 

= sup |e>v'^^(^--)(Ar^Af, + lV35-n|_gAvd^T^35-i 
> sup Af - i 

0<s<T ^ 


the bound 


sup {A 

'-0<s<T 


}>N 


-6 


<Po 


4’\Z) - J^’\z) > i 


< 2Te 


A^/log(Af)T ^K. 

Ni^ 


< 


2TC„ 

]\IK-1+3S ■ 


For any given /3 > 0, we can choose k := /3 -I- 1 — 35 so that the result 
takes the form 


sup {A 

0<s<T 


}>N 


-5 


< 


TC 

w 


(79) 


with C = 2C'„, as announced. 


□ 
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9 Controlling the mean field dynamics 

The previous theorem contains our main approximation result for the 
mean field dynamics. As explained in Section!^ two more steps remain 
in order to complete the proof of Theorem 14.41 First, we have to show 
that the solutions of the regularized Vlasov-Poisson equation (nni 
converge to a solution of the proper Vlasov-Poisson equation as the cut¬ 
off is lifted with V —>■ oo. Second we have to prove the approximation of 
the continuous Vlasov-density by the discretized version [^t,o{Z)]. 
The proof of ft follows the method of Loeper [13] . 

Proposition 9.1. Let /o satisfy the assumptions of Theorem \4.4\ Let 
ff^ and ft be the solution of the regularized, respectively the unregular¬ 
ized Vlasov-Poisson equation with initial datum /o. Then we have for 
p € [1, oo] and N > 3; 

Wp{ft^, ft) < N-^ (80) 

where Co depends on supt_^{||pf^||oo, ||Pt ||oo}- 

Proof. Let := p[ff] and p{ := p[ft] denote the charge density 
induced by f^ and ft, respectively. Let (p^ = {Q^,Pl^) be the char¬ 
acteristic flow of ff^. For the (unregularized) Vlasov-Poisson equation, 
the corresponding vector-field is not Lipschitz. However, as we assume 
the existence of a solution ft with bounded density pt, the mean field 
force k * Pt does satisfy a Log-Lip bound of the form \k * pt(x) — k * 
Pt{y)\ < Clx-2/1(1-Klog"(lx-J/D), where log"(x) = max{0, - log(x)}. 
This is sufficient to ensure the existence of a characteristic flow ipt s = 
{Qis^Pps) such that ft = 

Since Wp < Wq for p < q, it suffices to prove the statement for the 
infinite Wasserstein distance Woo- We consider 7ro(x, y) := fo{x)S{x — 
y) e n(/o, /o), which is the optimal coupling yielding Wooifi^, /t)|t=o = 
WooifoJo) = 0. Let TTt = ((/7(^o,V’t^o)#7ro. Then tt* e IV{ff^,ft)'it G 
[0, T) and we can consider 

D{t) : = TTt-esssup { Vlog(V)|x^ - y^ \ + W - y‘^\] 

= TTo-esssup {\/log(V)I( x) - Q{o(2/)I + \Pt%x) - P/o(2/)l} 
= /o-esssup {\/log(V)|g(^o(x) - QIo{x)\ + \Pt%x) - Ploi^)]}, 

which is an upper bound on Woo{ft^, ft) for N > 3. We compute: 
d^Dft) < /o-esssup|^/log N\pI^{x) - P^ (x)| 

+ I * Pt{Qfi.x)) -k*p{ (Q{ (x)) 11 

< y/log{N)D{t) 

-h/o-esssup Ifcf * pf (Qf (x)) -k* p{{Q{{x))\. 
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Now note that 


k * p{{Q{{x}) = [ k{Q{(x) - q) dp{{q) = [ k{Q((x) - q) dft{q,p) 

= [ HQ{(x) - Q{{z))dfo{z) 
jR3xR3 

and analogously for k^ * pf {Qf {x)). Hence, we can write the inter¬ 
action term as 

\ks *pf {Qt{x)) -k*p( {Q{{x))\ 

= f [k^{Q^^ix)-Q^{z))-k{Q{{x)-Q{{z))]dfo{z) 

jR3xR3 

<[ \k^iQi^{x)-Q^^{z))-k^{Q{{x)-Q{{z))\dfoiz) (82) 

jR3xR3 

+ f \k^{Q{{x) - Q{{z)) - k{Q{{y) - Q{{z))\dfo{z) (83) 

jR3xR3 

The second term (1551) can be bounded as 

J\k^{Qliy) - Q{(z)) - k{Q{{y) - Q{{z))\dfo{z) 

= j\ksiQ{iy) - 9) - k(Q{(y) - q)\p{{q)dq 

< ll*f - k\\i\\p{\\oo < dTrN-^\\p{\\^, (84) 

where we used the fact that k and k^ differ only on a set of radius 
N~^, so that ||fc^ — fc||i < ll^||g|<^-i111 = 47r7V“'^. The term (1551) can 
be bounded by hrst using the mean-value theorem as 

- Qt{z)) - kf (Q{{y) - Q{{z))\ 

< (iVfcf (Qf (x) - Qf (^))| + |Vfcf (Q{(x) - Q{ (0))|) 

•(lof (a;) - Q{{x)\ + \Q^{z) - Q{{z)\y 

Taking the integral with respect to dfo{z), we estimate \ Q^{z) — Q{{z)\ 
by the esssup and ||V/c^||i with Lemma [6.II This yields: 

dMl) < C'olog(A^)(|Qf (x) - Qf(x)| -k/o-esssup|Qf( 2 :)-Qf( 0 )|). 

with Cq := 2CiCp as before and Cp the uniform bound on the charge 
density from (1551) . Taking also the /q— esssup over the x-variable, we 
get the bound 

2Co \og{N) /o-esssup \Q^iz) - Q{ (z) | 

< 2CoV^ogiN)D{t). 

Putting everything together, we have found 

d^D{t) < 2Co\/log{N)D{t) -k 47r||pt||ooA^“^- 
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Using GronwalFs lemma and the fact that D{0) = 0, we get 


< D{t) < e‘2CoVi°g(^)iV-^ 

from which the desired proposition follows. 


□ 


A standard Gronwall argument yields the following result. 

Proposition 9.2. Let ip^ = {Q{t, ■), P{t, ■)) the characteristic flow 
of ffl defined by (1121) and the lift to the N-particle phase-space 

defined in m- Then we have for all p G [1, oo) and N > 3; 


Proof. For Z G let tto{x, y) G n(/r^, /o) and define tt^ = {pfi^ pfi)ffTTo G 
[^t,o{Z)], fl^). Note that both measures are now transported 
with the same flow. Set 


D 


ft) := J (^yiogfV|a;^ - y^l + d7r4(a;,y) 

RSxR® 

j (^y/log N\Qt{x) - Qtiy)\ + \Pt{x) - Ptiy)\y dTro{x,y) 


RSxR® 


We compute: 

^DP{t) = p j dTTo{x,y) (^y/logN\Qt{x) - Qt{y)\ + \Pt{x) - Pfiy)^ 
\Pfix) - Pt{y)\ + |fcf * pf (Qt(cr)) - fcf * pf'{Qt{y))\) 
Using again the Lipschitz bound 

\k^ * pf {Qfix)) - *p^{Qfiy))\ < CMN)\Qt{x) - Qfiy)\ 

for N > 3, i.e. log(A^) > 1, we arrive at the estimate 


D{t) < D{0) + Co\/\ogN J D{s)ds 


( 86 ) 


Hence by Gronwall’s inequality: 

Wp(p^[$t.o(Z)],/f) = < D{t) < i^(0)e*^'>^^. 

Taking on the right-hand side the infimum over all tto{x, y) G n(/i^, /o), 
Wp{p^[d>t,o{Z)],fr) < 


□ 

In view of (1851) . it remains to establish an upper bound on the typ¬ 
ical rate of convergence for Wp{pQ [Z], fo) —?> 0. (Note that, other 
than that, the result of Proposition liOl is actually deterministic.) For¬ 
tunately, we can rely for this purpose on recent, partcularly strong 
concentration estimates obtained by Fournier and Guillin, 2014 [3]. 
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Theorem 9.3 (Fournier and Guillin). Let f be a probability measure 
on R” such that 3k > 2p: 

Mkif) ■■= f \z\'"df{z) < + 00 . 


Let be a sample of independent variables, distributed ac- 

N 

cording to the law f and consider [Z] := ^ Szi ■ Then, for any 

i—l 

€ > 0 there exist constants c, C depending only on k, Mk{f) and e such 
that for all N >1 and ^ > 0; 


J)>i < CNiN^)-— + Cl 5 <i a{N, C), 


with 

{ exp(—cTV^^) ifp > n /2 

exp(-cjV( j^^4i/c) )^) ^fP = ^2 (87) 

exp(-ciVC”/^) */pe[l,n/ 2 ). 

With these large deviation estimates, we get the following. 

Corollary 9.4. Let p G [1, oo), 7 < min{5, i, and N > 3. Then 
there exists constants c,C > 0 such that 


3t G [0,r] : [^tAZ)]Jt) > (1 + 

( 88 ) 




where we use the notation 6 V 2p := max{ 6 , 2p}. 


Proof. By assumption in Theorem 14.41 there exists k > 2p such that 
Mkifo) < + 00 . Applying Theorem 1131 with ^ = N-P^,e = 
and the finite-moment condition ( 1 ), we get constants C, c > 0 such 
that 


Wp(/r^[Z],/o) >iV-^ <C(e 


_C^l-(6V2p)7 


-iW 


0 - 


Thus, with Proposition 19.21 we conclude 


3t G [0,r] : 


N 


< c{^ 






Adding the bound Wp{fA ft) < N ' 5 g* 2 CoVToiTv Proposition 19.11 

the statement follows. □ 


Now we have everything in place to complete the proof of Theorem l4.4l 
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Proof of Theorem 14.41 Let p G [1, oo), 7 < g and N > 3. We split 
the approximation into 


+ Wp{p^[<i>tAz)]Jt). 


From Corollarv l9.4l we get constants c,Ci >0 such that 

Fo\3tG[0,T]--Wp{p^[<i>tAZ)]Jt) > (1 + 


<cA 






In Theorem 14.21 we choose /3 = ^ ~ 1 and get a constant C' so that, 
together with Lemma [5111 


3t G [o,r]: wAAA^tAZAAA^tAZ)]) > n 


T-S 


< TC'N^ 


(89) 

, CqT + 1 '2 

for any N > ' . Putting both estimates together and choosing 

7 < min{i, J}, we have found that 


3t€[0A]:WAf^nZ]Jt)> 


N-^ + (1 + v'logiV)iV-T' e*2CoVi°g(^) 


(90) 


with C 2 '■= Cl + C". We can simplify this result by noting that 

gAfyiog(JV) < jyi-3S foj. ^ We shall thus demand N ^_ 

W := which yields > max{e, 2(l+fyIdglV) e 2 CoTViog(iV)| 

and conclude that 


»o[3t G [0,r] : Wp(/xf [Z],fy) > iV-^+i-3^ 




(91) 

□ 


10 Weaker singularities, open questions 


While the present paper focuses on the Vlasov-Poisson equation, the 
method presented here can also be applied to interactions with milder 
singularities, see [T]. For better comparison with other approaches, in 
particular the reference paper of Hauray and Jabin, 2013, |6], we shall 
state here the corresponding results without further proof. General¬ 
ization to higher dimensions would be straight-forward, as well. 

We use the characterization of force kernels introduced in Def. nn 

Theorem 10.1. Let a < 2. Let k satisfy a S°‘ condition and 
satisfy a Sf condition with the additional assumption ([9]) and 


6 < 


1 

1 -I- a 


(92) 
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Assume (for simplicity) that foGL^d x normalized 

to f fo = 1, has compact support and let f be the unique solution of 
the Vlasov equation with force kernel k. For Z G let H^[Z] the 

unique weak solution of the (regularized) Vlasov equation with force kf 
and initial data fj,Q [Z]. Then we have molecular chaos in the following 
sense: For any ^8 > 0,7 < inin{i,5} and T > 0 there exists constant 
C\^Ci such that 

Po[3t e [0,r] : Wl{^l^[Z]Jt) > N-^] < +TC 2 N-P. 

(93) 

This can be compared to the results in IH], where a statement similar 
to ((Ml) is derived for the case 1 < a < 2 with a cut-off of order 

5 < i minj(94) 

For a G [1, 2), the upper bound on S given by (1941) ranges between | 
and while our upper bound from (1921) ranges between ^ and In 
particular, it is interesting to note that the cut-off required in [5] is 
smaller than ours for a < | but larger for | < a < 2. This suggests 
that the probabilistic estimates presented here fare better for strong 
singularities - in the sense of admitting a significantly smaller cut-off 
- while the method proposed in [6] provides better controls for mild 
singularities. 

Most notably, Hauray and Jabin are able to treat the case 0 < a < 1 
with no cut-off at all by providing an explicit control on the minimal 
particle distance (in (p, g)-space, strictly speaking, while integrating 
the forces over small time-intervals). As it stands, our method re¬ 
quires in any case a regularization of the microscopic dynamics. Since 
it proves very effective in this setting, it would be interesting to inves¬ 
tigate if it can be extended - or possibly combined with the approach 
of [S] - to further reduce the cut-off or, ideally, dispense with it alto¬ 
gether for sufficiently mild singularities. 

Acknowledgements 

We thank Maxime Hauray, Michael Kiessling, Detlef Diirr, Martin 
Kolb, Ana Canizares, Samir Salem and Young-Pil Choi for valuable 
comments on earlier versions of the manuscript. 


29 



References 


[1] N. Boers and P. Pickl. On mean field limits for dynamical systems. 
Journal of Statistical Physics, pages 1-16, 2015. 

[2] W. Braun and K. Hepp. The Vlasov dynamics and its fluctuations 
in the 1 /N limit of interacting classical particles. Communications 
in Mathematical Physics, 56(2):101-113, 1977. 

[3] R. L. Dobrushin. Vlasov equations. Functional Analysis and Its 
Applications, 13(2):115-123, 1979. 

[4] N. Fournier and A. Guillin. On the rate of convergence in wasser- 
stein distance of the empirical measure. Probability Theory and 
Related Fields, 162:1-32, 2014. 

[5] F.A. Griinbaum. Propagation of chaos for the Boltzmann equa¬ 
tion. Archive for Rational Mechanics and Analysis, 42:323-345, 
1971. 

[6] M. Hauray and P.-E. Jabin. Particles approximations of Vlasov 
equations with singular forces : Propagation of chaos. To appear 
in Annales scientifiques de TEcole Normale Superieure, 2013. 

[7] E. Horst. Global strong solutions of Vlasov’s equation - neces¬ 
sary and sufficient conditions for their existence. Banach Center 
Publications, 19(1):143-153, 1987. 

[8] E. Horst. On the asymptotic growth of the solutions of the Vlasov- 
Poisson system. Mathematical Methods in the Applied Sciences, 
16:75-85, 1993. 

[9] M. Kac. Foundations of kinetic theory. In Proceedings of the Third 
Berkeley Symposium on Mathematical Statistics and Probability, 
1954-1955, volume vol HI, pages 171-197. University of California 
Press, 1956. 

[10] M. K.-H. Kiessling. The microscopic foundations of Vlasov theory 
for jellium-like Newtonian N-body systems. Journal of Statistical 
Physics, 155(6):1299-1328, 2014. 

[11] D. Lazarovici. The Vlasov-Poisson dynamics as the mean-field 
limit of rigid charges. Preprint: arXiv:1502.07047, 2015. 

[12] P.-L. Lions and B. Perthame. Propagation of moments and reg¬ 
ularity for the 3-dimensional Vlasov-Poisson system. Inventiones 
Mathematicae, 105:415-430, 1991. 

[13] G. Loeper. Uniqueness of the solution to the Vlasov-Poisson sys¬ 
tem with bounded density. Journal de Mathematiques Pures et 
Appliquees, 86:68-79, 2006. 

[14] H. Neunzert. An introduction to the nonlinear Boltzmann-Vlasov 
equation. In C. Cercignani, editor. Kinetic Theories and the Boltz¬ 
mann Equation, volume 1048 of Lecture Notes in Mathematics, 
pages 60-110. Springer, Berlin, Heidelberg, 1984. 

[15] H. Neunzert and J. Wick. Die Approximation der Losung 
von Integro-Differentialgleichungen durch endliche Punktmengen. 


30 



In R. Ansorge and W. Tornig, editors, Numerische Behand- 
lung nichtlinearer Integrodifferential - und Differentialgleichun- 
gen, volume 395 of Lecture Notes in Mathematics, pages 275-290. 
Springer, Berlin, Heidelberg, 1974. 

[16] K. Pfaffelmoser. Global classical solutions of the Vlasov-Poisson 
system in three dimensions for general initial data. Journal of 
Differential Equations, 95(2):281-303, 1990. 

[17] J. Schaeffer. Global existence of smooth solutions to the Vlasov- 
Poisson system in three dimensions. Communications in Partial 
Differential Equations, 16(8-9):1313“1335, 1991. 

[18] H. Spohn. Dynamics of Charged Particles and their Radiation 
Field. Cambridge University Press, Cambridge, 2004. 

[19] A.-S. Sznitman. Topics in propagation of chaos. In Ecole d’Ete de 
Probabilites de Saint-Flour XIX - 1989, volume 1464 of Lecture 
Notes in Mathematics, pages 165-251. Springer, Berlin, 1991. 

[20] C. Villani. Optimal Transport Old and New, volume 338 of 
Grundlehren der mathematischen Wissenschaften. Springer, 
Berlin, 2009. 


31 



